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Motivated by recent experiments on semi-dilute wormlike micelles, we study the early stages 
of the shear banding instability using the non-local Johnson-Segalman model with a "two-fluid" 
coupling to concentration. We perform a linear stability analysis for coupled fluctuations in shear 
rate j, micellar strain W_ and concentration cj> about an initially homogeneous state. This resembles 
the Cahn-Hilliard (CH) analysis of fluid-fluid demixing, though we discuss important differences. 
First assuming the homogeneous state to lie on the intrinsic constitutive curve, we calculate the 
"spinodal" onset of instability in sweeps along this constitutive curve. We then consider startup 
"quenches" into the unstable region. Here the instability in general occurs before the intrinsic 
constitutive curve can be attained so we analyse the fluctuations with respect to the time- dependent 
homogeneous startup flow, to find the selected length and time scales at which inhomogeneity first 
emerges. In the uncoupled limit, fluctuations in 7 and W_ are independent of those in <j>, and are 
unstable when the intrinsic constitutive curve has negative slope; but no length scale is selected. 
For the coupled case, this instability is enhanced at short length scales via feedback with <j> and 
a length scale is selected, consistently with the recent experiments. The unstable region is then 
broadened by an extent that increases with proximity to an underlying (zero-shear) CH demixing 
instability. Far from demixing, the broadening is slight and the instability is still dominated by 8^ 
and 8 W with only small 8(f). Close to demixing, instability sets in at very low shear rate, where it 
is demixing triggered by flow. 

PACS numbers: 47.50. +d Non-Newtonian fluid flows- 47.20.-k Hydrodynamic stability- 36.20.-r Macro- 
molecules and polymer molecules 



I. INTRODUCTION 

For many complex fluids, the intrinsic constitutive 
curve of shear stress E as a function of shear rate 7 is 
non-monotonic, admitting multiple values of shear rate 
at common stress. For example, Cates' model for semi- 
dilute wormlike micelles |jj predicts that the steady shear 
stress decrease above a critical 7 = 7ci- At very high 
shear rates, fast relaxation processes must eventually re- 
store an increasing stress 0]. See curve ACEG of 
Fig. 0. For the range j c i < 7 < j C 2 in which the stress 
decreases, steady homogeneous flow (Fig. ||a) is unsta- 
ble Q. For an applied shear rate 7 in this unstable range, 
Spenley, Cates and McLeish J3| proposed that the system 
separates into high and low shear rate bands (7h and 
je) with relative volume fractions satisfying the applied 
shear rate 7. (Fig. |]b.) The steady state flow curve 
then has the form ABFG. Within the banding regime, 
BF, a change in the applied shear rate adjusts the rela- 
tive fraction of the bands, while the stress E se i (common 
to both) remains constant. Several constitutive models 
augmented with intcrfacial gradient terms have captured 
this behaviour- |, |, g, g |t| . 

Experimentally, this scenario is well established for 
shear-thinning wormlike micelles (lj, [D| (li) . The steady 
state flow curve has a well defined, reproducible plateau 
E sc i. Coexistence of high and low viscosity bands has 
been observed by NMR spectroscopy JT^, 16, fl7| . 




Further evidence comes from small angle neutron scatter- 
ing p^_B^_g9 L_g0TJg l], and from flow birefringence 
(FB) p3|, \2A[ 25|, |26fl, which reveals a (quasi) nematic 



FIG. 1: Schematic flow curve for wormlike micelles. 



birefringence band coexisting with an isotropic one. The 
nematic band of FB has commonly been identified with 
the low viscosity band of NMR; but see [^7], |28| . 

Here we consider banding formation kinetics. Ex- 
perimentally ||, II, H, |g |l|, H ID, in rapid up- 
ward stress sweeps the shear rate initially follows the 
steady state flow curve (AB in Fig. |l|) before depart- 
ing for stresses E > E sc i along a mctastable branch 
(BC). When this branch starts to level off (hinting of 
an unstable branch for 7 ~ j c i) the shear rate finally 
"top-jumps" . Under shear startup in the metastable re- 
gion jg < 7 ~ 7 c i, the stress first rapidly attains the 
metastable branch BC (sometimes via oscillations), be- 
fore slowly decaying onto the steady-state plateau E se i 
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y= dv/dy 



FIG. 2: (a) Homogeneous shear rate and (b) banded profiles. 



via a "sigmoidal" envelope exp[— (£/tng) q ], with stretch- 
ing exponent a. The time scale tng greatly exceeds the 
Maxwell time r of linear rheology, but is smaller for shear 
rates 7 closer to 7 c i; for most systems the stretching ex- 
ponent is in the range 1.5 ~ a ~ 3. In the data of |33"| , 
for example, tng ~ (7 — ii)~ p an d a » 2 throughout the 
metastable range, with a crossover to a ~ 1 for 7 ~ Tel 
signifying onset of instability. In more dilute systems |31f| 
the onset of instability can be marked by a huge stress 
overshoot subsiding rapidly to E sc i via damped oscilla- 
tions. 

In the same experiments |HJ the stress overshoot coin- 
cided with pronounced concentration fluctuations that 
first emerged perpendicular to the shear compression 
axis, at a selected length scale 0(l/im). These fluc- 
tuations were attributed by the authors of Ref. ||l[] to 
the Helfand-Fredrickson [M coupling of concentration to 
flow. Although this mechanism has not, to date, been 
widely employed in the theory of shear banding - an im- 
portant exception being [ p5[ (discussed below) - it has 
been recognised elsewhere p4 p3C| . [37| [w| [Jq]. ffl, [II], 
fl3f as an important feature of two components systems 
with widely separated relaxation times (e.g. polymer and 
solvent). The slow component (polymer; for our pur- 
poses micelles) tends to migrate to regions of high stress; 
if the plateau modulus increases with polymer concen- 
tration (dG/d<fr > 0), positive feedback enhances concen- 
tration fluctuations perpendicular to the shear compres- 
sion axis (as seen in Ref. |3l[|) and (sometimes) shifts 
the spinodal of any nearby Cahn-Hilliard (CH) demixing 
instability [|(|. 

Further evidence for concentration coupling comes 
from the slight upward slope in the stress plateau BF 
of many wormlike micellar systems. This is most readily 
explained (in planar shear at least) by a concentration 
difference between the coexisting bands [Q, |3^] . Any cou- 
pling to concentration has important implications for the 
kinetics of macroscopic band formation, due to the large 
time scales involved in diffusion. 

In this paper we model the initial stages of banding 
formation in the unstable regime. We use the non-local 
Johnson-Segalman (d-JS) model || Q for the dynamics 
of the micellar stress, since this is the simplest tensorial 
model that admits a flow curve of negative slope. To 
incorporate concentration coupling, we use a two-fluid 
model |36l |45|, [l6| for the relative motion of the mi- 
celles and solvent. We perform a linear stability analysis 



(similar in spirit to the CH calculation for conventional 
liquid-liquid demixing) for coupled fluctuations in shear 
rate, micellar stress and concentration about an initially 
homogeneous shear state. We calculate the "spinodal" 
boundary of the region in which these fluctuations are 
unstable. We then consider startup "quenches" into the 
unstable region, predicting the length and time scales at 
which inhomogeneity first emerges. We also discuss the 
physical nature of the growing instability, according to 
whether its eigenvector is dominated by the flow vari- 
ables or by concentration. 

The two-fluid d-JS model to be studied in this paper 
was introduced and analysed briefly by us in a previ- 
ous letter |4q |. In this work we discuss more fully the 
model's origin and approximations, and give detailed nu- 
merical and analytical arguments supporting the results 
announced in Ref. E| . 

For simplicity we consider only fluctuation wavevec- 
tors in the velocity gradient direction, k = ky, and (sep- 
arately) the vorticity direction k = kz_. (The stability of 
the latter in fact turns out to be unaffected by shear in 
our model.) Indeed, the spinodal is commonly defined us- 
ing only these fluctuations for which k.x = (4^] (since 
any component in the flow gets advected onto the ve- 
locity gradient direction), though this restriction is less 
relevant inside the unstable region where fluctuations can 
grow on a time scale similar to that of advection. 

The paper is structured as follows. In section || we 
introduce the model and describe its intrinsic constitutive 



curves. In Sec. IV) we perform a linear stability analysis 
for slow shear rate sweeps along the intrinsic constitutive 
curve, to define the spinodal onset of instability. In Sec. |v| 
we consider shear startup "quenches" into the unstable 
region. We conclude in section VI. 



II. THE MODEL 

The existing literature contains several approaches for 
coupling concentration and flow [m| [37| |39| 
fll| . The two-fluid model considered by us follows closely 
that of Milner J|8| , although we extend his work slightly 
by including a Newtonian contribu tion to the micellar 
stress, for reasons discussed in Sec. IV C 1 . (Milner was 
mainly interested in slow shear phenomena, for which the 
Newtonian terms are unimportant.) 

The basic assumption of the two-fluid model is a sep- 
arate force balance for the micelles (velocity v m ) and 
the solvent (velocity v s ) within any element of solution. 
These are added to give the force balance for the centre 
of mass velocity 

v = <Pv m + (l-(f>)v a , (2.1) 
and subtracted for the relative velocity 

v ie i = v m -v s , (2.2) 
which in turn specifies the concentration fluctuations. 
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We give these dynamical equations in Sec. [IB below. 
First, we specify the free energy. 



A. Free energy 

In a sheared fluid, one cannot strictly define a free en- 
ergy because shear drives the system out of equilibrium. 
Nonetheless, for realistic experimental shear rates many 
of the internal degrees of freedom of a polymeric solution 
relax very quickly compared with the rate at which they 
are perturbed by the externally moving constraints. As- 
suming that such a separation of time-scales exists, one 
can effectively treat these fast variables as equilibrated. 
By integrating over them, one can define a free energy 
for a given fixed configuration of the slow variables. For 
our purposes, the relevant slow variables are the fluid 
momentum pv and micellar concentration cj) (which are 
both conserved and therefore truly slow in the hydrody- 
namic sense), and the micellar strain W (which is slow 
for all practical purposes). More precisely W_ is the local 
strain that would have to be reversed in order to relax 
the micellar stress: 



W = E.E 1 - 6, with <5r' = E.5r 



(2.3) 



where <5r' is the deformed vector corresponding to the 
undcformed vector Sr. 

The resulting free energy is assumed to comprise sep- 
arate kinetic, osmotic and elastic components: 



F = F K (v) + F°(4>) + F°(p£). 
The kinetic component is 



F K (v) = 



d 3 xpv 2 . 



(2.4) 



(2.5) 



The osmotic component is 



= / d 3 x [/(» + |(V^) 2 



dMi + ev)rw9)r, (2.6) 



where /" is the osmotic susceptibility and £ is the equi- 
librium correlation length for concentration fluctuations. 



The elastic component is 

F e {W) = \l d 3 xG((j))tr \K-log(£ + W)] (2.7) 
in which G(0) is the plateau modulus. 



B. Dynamical equations 

We now specify the dynamics. In any fluid element, 
the forces and stresses on the micelles are as follows: 
1. The viscoelastic stress a of the micellar backbone: 



g = 2(m + i).^=G(<t>)K. (2.8) 



2. The osmotic force 4>V_[SF° / '8<j)\ which induces con- 
ventional cooperative micellar diffusion and a "non- 
linear elastic force" <frV[5F c /8<f)}. 

3. A Newtonian stress 2<prj m rP from fast micellar 
relaxations such as Rouse mocfes, where 



D° =D - ±<STrL> , (2.9) 



and 



R m = I [Y« m + (V w m ) T ] . (2.10) 

We call r] m the "Rouse viscosity" (distinct from the 
zero shear viscosity of the total micellar stress), rj m 
is assumed to be independent of </), but is prefac- 
tored by the extensive factor ip. 

4. The drag force C(4>)llrei impeding the relative mo- 
tion of micelles and solvent. Scaling theory pcfl 
suggests C ~ 67r ?7^~ 2 where r\ is the mean viscosity 
defined in Eqn. 2.14 below. 



5. Stress due to gradients in the hydrostatic pressure 
P- 

The overall micellar force balance equation is thus: 



Pm^ (d t +H m .V) v m = V.G(c/>)m- + 2 V. HraR° m - CMErd - (2-11) 

Likewise, for the solvent we have the Newtonian viscous stress, the drag force and the hydrostatic pressure: 

Ps (l - 0) {d t + v s .V)v 8 = 2V. (I - <P)vsR° s + C(4>)v Tel - (I - 0)Vp. (2.12) 

These equations contain the basic assumption of "dynamical asymmetry" , i. e. that the viscoelastic stress acts only 
on the micelles and not on the solvent. Adding them, and assuming equal mass densities p m — p s — p, we obtain the 
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overall force balance equation for the centre of mass motion 



r 



in which I? and D° , defined analogously to -D 



r- — I = r- -i^= reX = m 

Eqns. [2.9| and [2.10] above, are the traceless symmetrised 
shear rate tensors for the centre of mass and relative ve- 
locity respectively. We have also defined 



and 



V = (f>r/m + (1 - ^)?7s 



rj = 4>(1 - 4>)(r] m ~ 77 s ). 



(2.14) 



(2.15) 



The equal an d op posite drag forces have cancelled each 
other in Eqn. 2.13 , which is essentially the Navier-Stokes' 
equation generalised to include osmotic stresses. The 
pressure p is fixed by incompressibility, 



0. 



(2.16) 



We attach a cautionary note to Eqn. 2.13 . The right 
hand side (RHS) is the net force acting on the fluid 



element. The LHS, however, equals the usual convec- 
tive inertial force plus the two terms in w rcl . We con- 
sider this to be an unsatisfactory aspect of the two-fluid 
model that is seldom acknowledged in the literature. One 
might argu e t hat t he separate advected derivatives of 
Eqns. [2.1l| and 2. IS should have v. 'Vv i in place of v i . 'Vv i 
(for i € m, s). This s till le aves the correction — v T ^v.V(t) 



on the LHS of Eqn. 2.13 and does not improve the ap- 
proximation. In this paper, however, we consider only 
small fluctuations about a homogeneous shear state (in 
which v rcl = 0), and the correction terms are truly neg- 
ligible. 



S ubtra cting the micellar and solvent Eqns. 2.11 



and 2.12 (with each predivided by its own volume frac- 
tion) gives the relative motion, which in turn specifies 
the concentration fluctuations: 



J 



{d t + v.V)(t>= -V.<Kl 



2 (1-^) 2 



A/? 1 

-V— + -V.G (<f>)W + 2Mi{D°) + 2A/" 2 CD° ,) 

Ash A\ ICl 



M\ and Mi are the following Newtonian terms: 

MUZ ) = ^Y.^° 

and 



rY.(l - <t>)D a 



^(D°J = ( J T 



1 



(2.17) 



(2.18) 



(2.19) 



We have neglected the inertia in 2; rol since it is small compared with the drag force Cl^roi [0- The osmotic contribution 
to the first term in the square brackets of Eqn. 2.17 shows that micelles diffuse down to gradients in the chemical 



potential, as in the usual CH description. The second term states that micelles diffuse in response to gradients in 
the viscoelastic stress. As described by Helfand and Fredrickon |34| , this provides a mechanism whereby micelles can 
diffuse up their own concentration gradient: the parts of an extended molecule that are in regions of lower viscosity 
will, during the process of relaxing to equilibrium, move more than those parts mired in a region of higher visosity 
and concentration. A relaxing molecule therefore on average moves towards the higher concentration region. This 
provides a mechanism whereby shear can enhance concentration fluctuations, and is the essential physics of the two 
fluid model. 

For the dynamics of the viscoelastic micellar backbone strain we use the phenomenological d-JS model jj| B : 



(d t + v m .V)W = a(D .W + W.D ) + (W.Q -fi W) + 2D 



W 



I 2 



v-it. 

= =m =m=' =m r (0) t((/))~ = 
I 



(2.20) 



where 2£ m = Vu m -(Vi; m ) T with (Vw m ) Q/3 = d a (v m )p. 
t(4>) is the Maxwell time. The length I could, for exam- 



ple, be set by the mesh-size or by the equilibrium cor- 
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relation length for concentration fluctuations. Here we 
assume the former, since the dynamics of the micellar 
conformation are more likely to depend on gradients in 
molecular conformation than in concentration. The equi- 
librium correlation length £ of course still enters our anal- 
ysis, through the concentration free energy of Eqn. 'Lt . 
Together, I and £ set the length scale of any interfaces. 
In the absence of concentration coupling this non-local 
term involving I is needed to reproduce a steady banded 
state with a uniquely selected stress ||, although other 
treatments [32| have used alternative forms for non- 
local terms that also give a uniquely selected stress. The 
slip parameter a measures the non-affinity of the molec- 
ular deformation, i.e. the fractional stretch of the poly- 
meric material with respect to that of the flow field. For 
\a\ < 1 (slip) the intrinsic constitutive curve in planar 
shear is c apable of the no n-m onoto nicity of Fig. H. Wc 
use Eqns. |2.I3| , |2.16| , |2. 17| and |2.20| as our model for the 
remainder of the paper. 



C. Flow geometry. Boundary conditions 



Parameter 


Symbol Q 


Value at <p = 0.11 


d log Q 
d log (j) 


Rheometer gap 


L 


0.15 mm 





Maxwell time 


T 


0.17s 


1.1 


Plateau modulus 


G 


232 Pa 


2.2 


Density 


P 


10 3 kgm -3 





Solvent viscosity 




l(r 3 kgm _1 s -1 





Rouse viscosity 




0.4kgm _1 s _1 





Mesh size 


I 


2.6xl0~ 8 m 


-0.73 


Diffusion coefficient 


D 


3.5 x lO^'m's" 1 


0.77 


Drag coefficient 


c 


2.4 x 10 12 kgm" 3 s~ 1 


1.54 


Correlation length 




6.0 x 10~ 7 m 


-0.77 


Slip parameter 


a 


0.92 






TABLE I: Experimental values of the model's parameters at 
volume fraction cj> = 0.11 (column 3). Scaling laws for the de- 
pendence of each parameter upon <f> (column 4). In most cal- 
culations we use the reference values of column 3 at (j> = 0.11, 
then tune <f> using the scaling laws of column 4. Only where 
stated do we allow the parameters to vary independently. 



We consider idealised planar shear bounded by infinite 
plates at y = {0, L} with (w,Vt/,VAi)) in the (x,y,z) 
directions. The boundary conditions at the plates are as 
follows. For the velocity we assume there is no slip. For 
the concentration we assume 



dy<f> = 8l<t> = 0. 



(2.21) 



which ensures (in zero shear at least) zero flux of concen- 
tration at the boundaries. Following ||, for the micellar 
strain we assume 



d y W aP = 0Va,/3. 



(2.22) 



Conditions 2.21 and 2.22| together ensure zero concentra- 
tion flux at the boundary even in shear. For controlled 
shear rate conditions (assumed throughout) 



dyj(y) = constant. (2.23) 



D. Model parameters 

Values for the model parameters are taken as fol- 
lows. We assume the solvent viscosity i] s and density 
p to be those of water. We take the plateau modulus 
G and the Maxwell time r from linear rheology (53) at 

= 0.11 on CTAB(0.3M)/NaNO 3 (1.79M)/H 2 O. We es- 
timate the Rouse viscosity r\ m from the (limited data 
on the) high shear branch of the flow curve of a closely 
related system [ p3| . The mesh size is estimated to be 

1 Ps (fcsT/G) 1 / 3 ]5C[| . In fact this form is only truly valid 
for a good solvent, although in the interests of simplicity 
we assume it to be good approximation even for systems 
closer to demixing. We take the diffusion coefficient D 



and the equilibrium correlation length £ from dynamic 
light scattering (DLS) data @ on CTAB/KBr/H 2 0, at 
a comparable micellar volume fraction. We calculate the 
drag coefficient £ = 6irr/^ 2 [Q. We fix the slip pa- 
rameter a — 0.92 by comparing our intrinsic constitutive 
curve in the semi-dilute regime to that of Cates' model 
for wormlike micelles {|]. We then have realistic values 
for all parameters, at cf> = 0.11 (table |). 

Exploring this large parameter space is a daunting 
prospect so we shall not, in general, vary the parame- 
ters independently. Instead we simply tune the single 
parameter cj>, relying on known semi-dilute scaling laws 
for the dependence of the other parameters upon (col- 
umn 4 of table ||) . For simplicity we assume that the slip 
parameter a is independent of cj>. We rescale stress, time 
and length so that G(c/> = 0.11) = 1, r(<f> = 0.11) = 1, 
and L = 1, where L is the rheometer gap (0.15mm) used 
in Ref. [ p3| . We also often eliminate p in favour of the 
Reynolds time Td = pL 2 /rj s . In total the model has 8 
scaled parameters. 



E. 



intrinsic constitutive curves 



In planar shear, the stationary homogeneous solut 
to Eqns. 2.13| , 2.1(\ , 2.17 and p. 20 for given 7 and 



10ns 
are 
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v lcl = and 



W xy = 



-yr(0) 



W vv = - 



1 + 6 7 2 t 2 (0)' 
1 bj 2 



where 



(1 + a) I + 67 2 
1 + a 

w zz = W xz = W yz = 0, 



1 



(2.24) 



(2.25) 



The total shear stress is the sum of the micellar stress 
and a Newtonian component: 



E(7,0) = G{4>)W xv + r ] {4>)^. 



(2.26) 



This defines a set of intrinsic constitutive curves £(7, <p) 
(dashed lines in Fig. ^). The criterion for the non- 
monotonicity of W xy to dominate the Newtonian term 
77(^1)7 and cause non-monotonicity in the overall stress 
£ is 77(0) < |G(</>)t (</>). As is reduced, therefore, the 
region of negative slope narrows, terminating in a "criti- 
cal" point at 4> = (f> c w 0.015. The same qualitative trend 
has been seen in CPCl/NaSal/brine 12 . 



o Infinite C, 

n D=2.6e-4 

« D=2.6e-5 

a D=2.6e-6 



10 

xy .--'^ 
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FIG. 3: intrinsic constitutive curves for (j> = 
0.11,0.091,0.072,0.053,0.034,0.015 (dashed lines, down- 
wards). Spinodals for: the uncoupled limit £ — > 00 (o); 
coupled model with D((j> — 0.11) taken from the DLS data 
(table |) (□); coupled model with D artificially reduced 
(0,A). Inset: enlargement at high 7. 



III. UNCOUPLED LIMIT; INSTABILITIES; 
POSITIVE FEEDBACK 



In the limit of infinite drag, i.e. £ — > 00 at fixed 
/"(</>), the relative motion between micelles and solvent 
is switched off, disabling concentration fluctuations. In 



the slightly different limit of £ 
diffusion coefficient 



00 at fixed micellar 



D = 



\2 fll 



c 



(3.1) 



2.171 then reduces to the 
Indepen- 



the concentration still fluctuat es, b ut independently of 
the rheological variables. Eqn. 
CH equation (with a ^-dependent mobility), 
dently of Scf>, the shear rate and micellar stress together 
obey uniform-0 d-JS dynamics ||, |5^, |56|. Accordingly, 
two separate instabilities are possible: 

1. Demixing instability 

For D < 0, concentration has its own CH demixing 
instability, gove rne d primarily by the free energy 
defined in Eqn. 2.6. In this work we consider only 



flow induced instabilities, for which D > 0. 

2. Mechanical instability 

For shear rates where the intrinsic constitutive 
curve has negative slope d'E/d'j < 0, fluctuations 
in shear rate and micellar stress have their own 
shear banding instability, which for convenience we 
call "mechanical" . However, we emphasize that the 
term "mechanical instability" is often misused in 
the literature. For example, the origin of banding 
in semi-dilute micelles is usually described as purely 
"mechanical" . In contrast banding in concentrated 
systems, nearer an equilibrium isotropic-ncmatic 
(I-N) transition, is generally described as a flow 
induced perturbation to this I-N transition. How- 
ever, such a transition possesses an essentially iden- 
tical "mechanical" instability so there is no sharp 
distinction between these scenarios. Our approach 
is more obviously relevant to semi-dilute systems 
since we rely purely on non-linear coupling between 
shear and W_ to induce instability, rather than on 
a perturbation of a non-linear free energy F C (W). 
Accordingly, but for convenient nomenclature only, 
we refer to this instability as "mechanical" . 

For finite drag, these instabilities mix. Consider, for 
exam ple, the second term in the square brackets of 
Eqn. 2.17, which encourages micelles to move up gra- 
dients in the micellar backbone stress. If the micellar 
stress increases with concentration (dG/d(f> > 0), positive 
feedback occurs leading to shear enhancement of concen- 
tration fluctuations J34|, [}?], [38| |3^| or (from the op- 
posite extreme) concentration-coupled enhancement of a 
flow instability. Instability can then occur even if D > 
and d'E/d A / > 0, so the domain of (what was the purely) 
mechanical instability is broadened relative to that in 
which d'E/d'j < 0. For systems far from a zero-shear CH 
demixing instability (i.e. for D 3> 0) the mixed instabil- 
ity is "mainly mechanical" (dominated by S'y, 6W). For 
systems that in zero shear are already close to demixing 
(D ~ 0) the instability sets in at low shear rate, where it 
is essentially demixing triggered by flow (dominated by 
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5(f)). Of course, as in the above discussion of mechani- 
cal vs.I-N instabilities, there is no sharp distinction be- 
tween these two extremes: the model captures a smooth 
crossover between them. 

Enhancement of flow instabilities by concentration 
coupling was first predicted by the remarkable insight of 
Schmitt et al. |}5|. The feedback mechanism described 
in the previous paragraph for our model corresponds to 
the direct assumption by Schmitt et al. [ fj5| of a chemical 
potential \x — /i(7). However this is only truly equivalent 
to our approach if the viscoelastic stress ^(7) can ad- 
just adiabatically (assumed in |3j|). Wc find below that 
the dynamics inside the spinodal are dictated by the rate 
of micellar stress response. The spinodal is unaffected, 
since response here is adiabatic by definition. Schmitt et 
al. also predicted an instability for negative feedback, 
but concluded it to be similar in character to a pure 
mechanical instability in which concentration plays no 
role. In our model, negative feedback would correspond 
to dG/dcj) < 0; here we consider only positive feedback 
(in the language of Ref. ||, C > 0). 

In this work, we model the onset of instability for two 
different flow histories. The first (section [V) assumes 
an initial state on the intrinsic constitutive curve, and is 
used to define the "spinodal" limit of stability in sweeps 
along this flow curve. This is analogous to defining the 
spinodal of a van der Waals fluid via quasistatic compres- 
sion, and suffers the same practical ambiguity that finite 
fluctuations can cause separation/banding via metastable 
kinetics before the spinodal is reached. The second his- 
tory (section |v|) considers a startup "quench" into the 
unstable region, and is (essentially) the counterpart of a 
temperature quench into the demixing regime of a van der 
Waals fluid. The analysis here is complicated by the fact 
that the fluctuations emerge against the time-dependent 
startup flow. 



IV. INITIAL CONDITION ON INTRINSIC 
CONSTITUTIVE CURVE; SLOW SHEAR RATE 
SWEEPS 

A. Linear analysis 

We encode the system's state as follows, 



u(r_,t) is real. Linearising Eqns. 2.13, 2.16, 2.17 and 2.2C 
in these fluctuations, we find 



(4.1) 



in which the e are dimensionless unit vectors |57j. Con- 
sider fluctuations about a mean initial state u that is on 
the intrinsic constitutive curve: 



u(r, t) = u + 3lLk(t) ex P (ih-L) 



(4.2) 



where u = [7, W, <t>\ with W(j) the stationary homoge- 
neous solution given by Eqn. p.24| . The sum in Eqn. 4.2 
covers positive and negative fc, with u_ k = u k ensuring 



d t Su k (t) =K k -Su k (t) 



(4.3) 



The stability matrix M fc depends on the model param- 
eters and the initial state u = [7, W \ <f>]. Its eigenmodes 
are determined by 



(4.4) 



where a is the mode index. The functions cuk,a versus k 
define a multi-branched dispersion relation. For the ini- 
tial state u to be stable, all dispersion branches must be 
negative. A positive eigenvalue u>k, a indicates an unsta- 
ble mode that grows exponentially in time with relative 
order-parameter amplitudes specified by the correspond- 
ing eigenvector v fe a . In an upward sweep along the in- 
trinsic constitutive curve, the lower spinodal lies where 
the eigenvalue ujk' with the largest real part (maximised 
over k and a) crosses the imaginary axis in the posi- 
tive direction. The upper spinodal is defined likewise, for 
sweeps towards the unstable region from above. Strictly, 
only harmonics of the gap-size are allowed; but in order 
to define the spinodal independently of rheometer geome- 
try we allow arbitrarily small wavevectors. For any shear 
rate between the spinodals, the dispersion relation is pos- 
itive for some range of wavevectors. Typically, we find 
just one unstable dispersion branch u! k (although we com- 
ment below on an exception, for some model parameters, 
at very high shear rates). 

Below we give numerical results for this unstable 
branch, focusing on any global maximum, which would 
indicate a selected length scale fc* _1 at which inho- 
mogeneity emerges most quickly. We also give results 
for the unstable eigenvector v fe . at this maximum. As 
noted above, we consider only fluctuations of wavevectors 
k = ky and (separately) k — kz. In fact the stability of 
fluctuations k = kz turns out to be unaffected by shear 
in our model (see section 1VC3) so we study in detail 
only k = ky. 

For k = ky, Sv y — by incompressibility. Fluctuations 
in the remaining variables decouple into three indepen- 
dent subspaces: 

• 61 = [ik8v x = 6-y, SW xy , 8W XX , SW yy ,8(f>\ 

• 6 2 = [ikSv z , W xz , W yz ] 

• 63 = w zz . 

In all unstable regimes, for this flow history, only Si 
is unstable. (In startup at high shear rates, S2 can go 
unstable; however it is always less unstable than 61 in 
the relevant time window; see section ^.) Accordingly, 
we focus on 61. For convenience, we change variables to 



1 + a w 
—2~ W vv 



(4.5) 
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and 



a — 1 1 + a 
Y = ——W xx —W v 



(4.6) 



for the instability in this subspace ©i. In some regimes 
we also give qualitative analytical results, obtained using 
the following simplified stability matrix of &%, 



In sections IV B and IV C we give our numerical results 



M fc = 





r\k 2 


k 2 





G'W xy k 2 


\ 


8"f 




VsTd 


VsTd 




VsTd 








1 + Z 


-l-l 2 k 2 


7 


W xv t' 




sw xy 




-bW xy 


—b'j 


-1 - l 2 k 2 


Zt' 




sz 


V 








k 2 /c 


~bk 2 {\ + ^ 2 k 2 ) 


) 


5<p 



(4.7) 



with 



and 



(l + o) 



D = D - 



za 



(4.8) 



(4.9) 



[Z is negative so D > D.) Matrix 4.7 is exact in the 
uncoupled limit ( — > oo. For finite £ it contains several 
approximations [g8| (most notably neglecting SY) and so 
underestimates the growth rate of the coupled instability; 
however the qualitative trends are unaffected. In some 
places below we further neglect terms of order 77. This 
is only valid for concentrations not to near the critical 
concentration </> c and shear rates not too far above the 
lower spinodal, so that 777 <c GW xy . In any case, startup 
at higher shear rates is too violent to study experimen- 
tally M. 



B. Results: uncoupled limit 

In the limit £ — > 00 at fixed fluctuations in the 

rheological variables decouple from those in concentra- 
tion and the stability matrix is exactly 



M , = Mm 



mechanical 



-Dk 2 {l + t 2 k 2 ) 



(4.10) 



in which M is the upper- left 3x3 "mechanical" sector 
of matrixTO. The three elements represented by the 
dash in matrix 4.10 are non-zero, but irrelevant since all 



off-diagonal elements in the bottom row are zero. In this 
section, we give results for the spinodal and dispersion 
relation in this uncoupled limit. 



1. Spinodal 



For each of a range of concentrations, we calculated 
spinodals numerically: see the circles in figures I The 
unstable region coincides with that of negative constitu- 
tive slope dEi/dj, as expected, and vanishes at a "critical 
point" 4> c w 0.015, as in the experiments of Ref. p2|| . 

A nalyt ically, the eigenvalues of the stability ma- 
trix 4.10 obey the quartic equation: 



bco 2 k 



■ CLU k 



(4.11) 



where d = Det M . The roots of any polynomial with 
real coefficients are either real, or complex-conjugate 
pairs. This gives two possibilities for the spinodal. First, 
the root with the largest real part could be zero, imply- 
ing the onset of a monotonically growing instability; this 
corresponds to the 3-sub-space (of the 4-space spanned 
by a, 6, c, d) for which d = 0, a > 0, b > 0, c > 0. Alter- 
natively, the root could be one of a pure imaginary pair, 
implying the onset of growing oscillations; this also cor- 
responds to a 3-sub-space though not (in general) defined 
simply by one of the a, 6, c, d axes. For the parameters 
considered, we have mostly found the first case.[^4| Ac- 
cordingly, our analysis hereafter considers only this first 
case for which the spinodal is given by d — Det M_ k = 0. 

(In all parameter regimes studied, this automatically en- 
sures a > 0, b > 0, c > 0.) It can further be shown that 
Det M < in the unstable region, i.e. 



DT) M > 



(4.12) 
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in which 



2) 



M 



Dct M 



M 



r\ s T d 



{- V (l + b^ 2 )-(l + Z) + bW xy ^} 



VsTd 



-(1 + b^ 2 ) 



(4.13) 



To a good approximation we have neglected the intcrfa- 
cial terms in calculating the spinodal: they merely cut 
off the dispersion relation at short length scales without 
affecting the sign of the maxi mum g rowth rate. The term 
in the square bracket of Eqn. 4.13 is always positive, so 



the condition for instability is finally just 

-D-~ > 0. 

07 



(4.14) 



For an increasing flow curve ^jp > 0, CH <j>— demixing 
can occur for D < 0. As noted above, in this paper we 
consider only D > 0, for which the unstable region is 
^? < 0, as shown numerically. Here, the instability oc- 



curs in the upper 3x3 subspace of the matrix 4.10| , and 
is purely mechanical. Although the normal stresses (en- 
coded by Z = Z(W XX ,W VV )) have apparently cancelled 
from Eqn. 4.14, they in fact play a crucial role, as follows. 



The origin of the instability is the term ■^rbW xy A f in the 



curly braces of Eqn. |4.13| . In this term, is the pref- 



actor to SW xy in the S'j equation, and states that a local 
increase in W xy causes a (diffusive) decrease in 7. The 
remaining factor feeds back positively, by stating that a 
decrease in 7 in turn tends to increase W xy , consistent 
with the negative slope in the constitutive curve. How- 
ever this factor itself comprises two subfactors, each of 
which describes a mechanism that explicitly involves the 
normal stress, Z. The first, —bW xy , states that the de- 
crease in 7 causes an increase in Z. The second, 7 states 
that this increase in Z causes an increase in 7, thereby 
completing the positive feedback. This role of normal 
stress was not considered in early studies of mechanical 
instability Note finally that the absolute values of 
the normal stresses are important, not just the difference 
W yy — W xx : the trace of the micellar contribution to the 
stress tensor is not arbitrary. 



2. Dispersion relation 

Before discussing the dispersion relation we make the 
following cautionary remark. This dispersion relation de- 
scribes fluctuations about a homogeneous state on the in- 
trinsic constitutive curve. While we correctly employed 
it to define the spinodal boundary of instability, it is less 
useful inside the unstable region since one cannot pre- 
pare an unstable initial state. Indeed, startup quenches 
into the unstable region in general go unstable long be- 
fore the intrinsic constitutive curve can be attained (see 



section |y|). However, the main features of this dispersion 
relation do still appear in their time dependent coun- 
terparts of startup flow. Our motivation for discussing 
them here is to gain early qualitative insight without the 
complication of time-dependence. 

For this uncoupled model (with this initial condition) 
we o nly ob serve one positive dispe rsion branch, shown in 
Figs. |(a]] for = 0.11 and |(a}| for = 0.02. Strictly, 
only harmonics k — mi of the gap size L = 1 are allowed. 
However we show k < it as well, because for some systems 
the features of this domain (disc ussed below) could lie in 
the allowed region k > n. Fig. [4(d)| contains the same 



data as 4(a), but enlarged on shear rates near the lower 
spinodal: this is the only regime in which startup kinetics 
have been studied experimentally since they become too 
violent at higher shear rates ^9| . 

For a given unstable applied shear rate 7, the growth 
rate lo tends to zero as k — > and as k — > 00, with a broad 
plateau in between. This can be understood via the fol- 
lowing analytical r esults obtained from the chara cteris tic 
equation of matrix 4.1C , and schematised in Fig. 3(a)| . 



(i) Reynolds regime k — > 0. Here we find 



k 2 
dj r] s T d 



(4.15) 



This is marked as a dashed line in Fig. 6(a), and 
agrees well with the numerical data. Here, the in- 
stability is limited by the Reynolds rate at which 
shear rate (conserved overall) diffuses a distance 
0(l/k): the micellar stress responds adiabatically 
in comparison. 

(ii) Non-conserved plateau regime. At these shorter 
length scales (but still with k 2 l 2 <C 1) the growth 
rate is instead limited by the Maxwell time on 
which the micellar backbone stress evolves (the 
Reynolds number is then effectively zero) . Because 
micellar stress is non conserved, the growth rate is 
independent of k: 



D 



M 



1 



l + Z {l + Z) 2 07 



■0(r],fj) 



with 



Dm = B 



M ' 



k 2 ' 



(4.16) 



(4.17) 



which is marked as a dashed line (also incor porating 
the interfacial regime, below) in Fig. |S(a) . 



(iii) Interfacial cutoff. The dispersion relation is cut off 
once kl — 0(1) by th e relu ctance to form interfaces. 
Here, w follows from [4.16] with uj — ► uj + Ik. 



The crossover between the first two regimes occurs at a 
length scale much greater than the interfacial cutoff, giv- 
ing a broad intermediate plateau. The maximum in ui(k) 
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log(k/ji) 2 

o 

(b) D(4> = 0.11) = 2.6 x 1CT 4 . 




log(k/jt) 2 

o 

(c) D(4> = 0.11) = 2.6 x 10" 6 . 




o 

log(k/it 




o 

log(k/Jt) 




o 

log(krt) 



(d) C -+ oo (detail). 



(e) D(<j> = 0.11) = 2.6 X 10" 4 (detail). (f) D(<f> = 0.11) = 2.6 X 10~ 6 (detail). 



FIG. 4: Positive (unstable) dispersion branch at <f> = 0.11. 4(a) and 4(d) are for the uncoupled model; 4(b) and l(e] are for the 
coupled model in which all parameters assume the experimental values of table |l| (spinodal given by Ds in Fig. \l(c] and p(f)| 
are for a coupled model in which D((f>) is artificially reduced (spinodal given by As in Fig. 0). For each vertical pair of graphs, 
the bottom is an enlargement of the top one, at shear rates near the lower spinodal. In each subfigure, the white space defines 
(7, k) values for which all dispersion branches are negative. 



is very shallow and its length scale exceeds the system 
size for the experimental systems considered here. There- 
fore fluctuations grow equally quickly at all length scales 
from the system size down to the interface width, and 
there is no selected length scale. In a previous work |5|] 
that considered the mechanical instability of a simple 
model (with no concentration coupling), a wavevector 
was apparently selected. However, there the viscoelas- 
tic stress was assumed to respond adiabatically so the 
intermediate plateau was absent. 



C. Results: coupled model 

For finite drag, fluctuations in the mechanical vari- 
ables are coupled to those in concentration S(f> via two 
main mechanisms. The first (already discussed briefly) 



involves the 2nd term in the square brackets of Eqn. 2.17, 
which decrees that concentration diffuses in response to 
gradients in W vv at rate tx 1/C; the elastic part of the 



stress (Eqn. 2.13) then increases in proportion to G'((j)) = 
dG((f>) /dtp, giving positive feedback oc G'/£. The second 
(neglected in our analytical work, as already noted ||58|) 
comes from the elastic contribution V(<5F C / 'S<fi) to the 
first term in the square brackets of Eqn. p,17| . 



1. Spinodal 

The mechanical instability is enhanced by this con- 
centration coupling. For the model parameter values of 
table I the spinodals are shifted only slightly (squares 
in Fig. H). However the shift increases near an under- 
lying CH demixing instability, as illustrated by reducing 
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log(k/Ji) 2 ^-^/io log(k/jt) 2 

4 o ' o 

(b) D(4> = 0.11) = 2.6 x 1CT 4 . (c) D(4> = 0.11) = 2.6 x 10" 6 . 

FIG. 5: Unstable dispersion branch at <f> = 0.02. 5(a) is for the uncoupled model; 5(b) is for the coupl ed m odel in which at 
4> = 0.11 all parameters assume the experimental values of table | (spinodal is given by Ds in Fig. |3); 5(c) is for a coupled 
model in which D((f> — 0.11) is artificially reduced (spinodal given by As in Fig. W). In each subfigure, the white space defines 
(7, k) values for which all dispersion branches are negative. 




FIG. 6: Illustration of the various dispersion regimes discussed in the text, (a) for the uncoupled mechanical instability (b) 
for the coupled model. (b,left) is for a shear rate that would be stable in the uncoupled limit £ — > 00; (b, right) is for a higher 
shear rat e that is in side the uncoupled mechanical spinodal. The dashed lines are t he ap proximate analytical asymptotes given 



in Eqns. 4.15, 4.16 



and 4.21 



The dashed arrows show the approximate k* of Eqn. 4.24. 



D(4> — 0.11) at fixed coupling G'/C (diamonds and tri- 
angles in Fig. |3|). This is intuitively obvious: when D 
finally goes negative (not shown) demixing must occur 
even in zero shear. In the opposite extreme D —> 00, 
the uncoupled limit is recovered (circles). The shifts in 
the lower spinodal have important implications for fast 
upward stress sweep experiments, since "top" jumping 
should in this case occur before the maximum of the un- 
derlying flow curve is reached. 



that qualitatively reproduces the shifts in the lower spin- 
odal (found by setting ojk = in the cha ract eristic equa- 
tion of the approximate stability matrix T7) is 



DD M + -Df > 

c 



(4.18) 



An approximate analytical condition for instability 



in which S)m is the mechanical determinant already de- 
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fined in Eqn. 4.13 and Dp is a "feedback determinant", 



S F = 





1 + Z 



-bW.„ 



■<:<) 



k 2 


G'W xy k 2 






-1 


W xy r' 


— 67 


Zt' 



0(77,77) 



-G'W X 



■G'W xy {-67(1 + Z)- bW xy ) + 0(77,77) 



7? s r s 



-(I + 67 2 ) 



^£ + 0(77,77), (4.19) 
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where dZ/dj < 0. (The terms in r' cancel each other.) 
As for the uncoupled model, the interfacial terms have 
been neglected in locating the spinodal. Our final condi- 
tion for instability is thus 



f ) d^ + G^W xy dZ <Q 
dj C dj 



(4.20) 



which reduces to the uncoupled condition 4.14 for £ — * 00 
at fixed D, as required. The size of the second term 
above (which encodes feedback) relative to the "diag- 
onal" product of uncoupled instabilities (first term) is 
set by G'/(D() ~ G'/f", i.e. the ratio of the "feedback 
elasticity" G" to the osmotic elasticity /". The kinetic 
coefficient £ has cancelled from this ratio, since th e in- 
stability occurs adiabatically at the spinodal. Eqn 4.2C 
corresponds to Eqn. 24 in the paper of Schmitt et al. [|35[. 

On the basis of these results, we classify systems into 
two basic types. 

• Type I systems are far from a CH demixing insta- 
bility (D ^> 0). The mechanical spinodal is shifted 
only slightly by concentration coupling. 

• Type II systems are close to a CH instability (D ~ 
0). The mechanical spinodal is strongly perturbed 
by concentration coupling. 

Correspondingly, we anticipate two types of instability 
(with a smooth crossover in between): 

• Type A instabilities, which are essentially mechan- 
ical (eigenvector mostly in Sj, SW) but perturbed 
by coupling to 5<fi. These are expected in all type I 
systems; and in type II systems for shear rates well 
above the lower spinodal. 

• Type B instabilities, which arc essentially CH in 
character (eigenvector dominated by S(f>), but per- 
turbed by a coupling to flow. These occur in type 
II systems at shear rates just inside the lower spin- 
odal: see Refs. SI M M M M Eol . 



This intuition is confirmed by the results given in sec- 
tion below. 

The results of Fig. || also reveal a second lobe of insta- 
bility that appears at high shear rates for small values of 



D. However its existence and location are highly sensi- 
tive to the choice of model parameters and to the precise 
details of model definition: it appears much more readily 
and extends to much higher shear rates if the Newtonian 
contribution to the micellar stress is not included. Its 
eigenvector is overwhelmingly dominated by 5^. It is as- 
sociated with two complex eigenvalues with equal positive 
real parts. We do not study this instability in detail, but 
return in section VI to discuss its potential implications. 



Z A a) C, infinite (uncoupled limit) 



^ A b) ^ finite (coupled model) 



FIG. 7: Sketch of the unstable (dashed) region of a mechanical 
instability a) decoupled from or b) coupled to concentration. 
As discussed in the text, concentration coupling broadens the 
region of instability, and can sometimes cause a new region of 
instability to develop in the high shear rate branch. 

The effect of concentration coupling in broadening the 
region of instability is schematized in Fig. [7], along with 
the possibility of a further region of instability in the high 
shear rate branch, discussed in the previous paragraph. 



2. Dispersion relation 

We now discuss the dispersion relation of the c oupled 
model (though the cautionary remark made in Sec. IV B 2 
above for the uncoupled model still applies). We focus 
mainly, and firstly, on shear rates above, but quite close 
to, the lower spinodal, since this is the regime that is 
studied experimentally. Comparing the dis persi on rela- 
tion for the pure mechanical instab ility (Fig. 4(a) ) to that 
for a coupled model of type I (Fig. 4(e)), we see that con- 
centration coupling enhances the mechanical instability at 
short wavelengths, thereby selecting a length scale. We 
discuss this length scale in more detail below. At long 
wavelengths the plateau of the uncoupled instability is 
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still apparent (provided dTi/d A f < 0) and unperturbed; 
for dT,/d A f > this plateau disappears to leave only the 
diffusive, concentration-coupled bump. The dispersion 
relation for a system closer to type II ( D red uced by a fac- 
tor 100, at fixed G'/C) is shown in Fig. [4(f)| : the enhance- 
ment at long length scales is much more pronounced, 
corresponding to the greater spinodal shift (triangles of 
Fig. |3|). However the mechanical plateau (present when 
d'S/d'y < 0) is still unperturbed at lon g len gth scales 
(though indiscernible on the scale of Fig. 4(f)). 

The overall dispersion shape (the same in type I and II 
systems) is captured by analysing the simplified stability 
matrix 



4.7. We consider two separate cases: 



(a) Shear rates above the lower spinodal of the coupled 
model but which a re sta ble in the uncoupled limit 
(d'S/d'y > 0; Fig. |(b| left). Here, we find the 
following regimes: 



(i) Diffusive regime fc — > in which 



(4.21) 



This is marked as a dashed line in Fig. 
|6(b)| (left), and underestimates the exact result 
because it neglects feedback via the elastic free 
energy F c pm . The growth rate in this regime 
is limited by the rate at which matter diffuses 
a distance 0(l/fc): momentum diffusion and 
micellar strain response are adiabatic in com- 
parison. Note that for larger shear rates for 
which dS/c?7 < (discussed in (b), below) 
Eqn. 4.21 is negative, so this branch is ab sent 
from the instability (compare Figs. |6(b)| left 
and right). 

(ii) Non-conserved "plateau" regime. For larger fc, 
the rate at which the non-conserved micel- 
lar strain can respond (even within concentra- 
tion enhanced dynamics) is the limiting factor; 
concentration diffusion becomes adiabatic in 
comparison. If the eventual interfacial cutoff 
in the dispersion relation once l 2 k 2 = 0(1) or 
£ 2 fc 2 = O(l) were absent we would then see a 



non-conserved fc- 
in which 



independent plateau regime 



M 



S) F 



l + Z 



with 



2>i = 2>i 



fc 2 ' 



G'bW 2 xy 
DC 



ie F,M. 



(4.22) 



(4.23) 



However, for the systems of interest to us 
the low-fc crossover to this regime is not well 
separated from the interfacial cutoff and the 



plateau is replaced by a rounded maximum at 
(k*,uj*) (thus defined) where This 
maximum selects a length scale fc* -1 . 

(iii) High fc interfacial cutoff. The dispersion rela- 
tion is cut off by interfaces once fc 2 / 2 = 0(1) 
or fc 2 £ 2 = 0(1). I and £ are of similar order 
for the systems of interest to us. 

An estimate for the selected wavevector fc* can be 
be obtained by expanding about lo sa ui p i to find 



, D(l + Z)l 2 G'bWlf ' 



(4.24) 



DC 



M 



M 



where 



Sm = S>m-^ p i(1 + ^) <0. 



(4.25) 



This fc* is marked by a dashed arrow in Fig. |6(b) , 
and agrees reasonably with the numerics. As in the 
conventional CH instability, fc* — * at the spinodal 
boundaries (where oj* — ► 0). This is not visible in 
figures [| and 0, because fc* only starts to diminish 
appreciably for indis cernibly small w* on our scale. 
Note that Eqn. |4.24 doesn't reproduce the selected 
wavevector of standard CH theory in zero shear, 
since phase separation is still affected by coupling 
of composition to viscoelastic effects [[KJ even in 
this limit. 

(b) For higher shear rates that would have been un- 
stable even in the uncoupled limit, the disper- 
sion relat ion develops a shoulder at small fc: see 
Fig. |6(b)| (right). As noted above, this is just the 
large length scale part of the pure mechanical dis- 
persion branch (Sec. IV B), comprising a Reynolds 



regime and a mechanical non- conse rved regime. 
(See regimes (i) and (ii) in Fig. 3(b) (right).) The 
growth rate here is much faster than diffusion so 
concentration is absent from the eigenvector. At 
shorter length scales, concentration can keep pace 
and is included. For shear rates that are not too 
deep inside the unstable region, the dispersion rela- 
tion then rises to the rounded plateau estimated by 
Eqn. 4.22 (regime (iii) of Fig. |6(b)| ( right)) before fi- 
nally being cut off by interfaces (regime (iv)). The 
maximum at fc* is again estimated by E qn. 4.24 
(marked by the dashed arrow in Fig. [6(b)l(rightyr 

The preceding analysis captures the qualitative fea- 
tures of the dispersion relations in many regimes. How- 
ever, some more exotic effects are apparent in Figs [4(b)| 
and 4(c) for shear rates well above the lower spinodal. 



For 20 ~ 7 ~ 80, concentration coupling gives negative 
feedback at short length scales. The origin of this (not 
included in our above analytical treatment) is that the 
velocity advecting the micellar backbone strain is not the 
centre of mass velocity v (as the above analytical work 
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assumed) but the micellar velocity v m = v + (1 — 4>)v 



A fluctuation 6W, 
and therefore in 



re] ' 



as 



yy 



in general causes a fluctuation in 
When included in the advective 



term, this feeds back negati vely on W yy . At still higher 
shear rates 7 > 80 in Fig. 4(c), the dispersion relation 
has a pronounced ridge corresponding to the high shear 
rate lobe discussed above and schematised by the right 
hand dashed line of Fig. 0b. 



3. Fluctuations in the vorticity direction 

In the uncoupled limit £ — > 00, the mechanical sub- 
space is stable with respect to vorticity fluctuations at 
all shear rates, while concentration has the usual CH 
demixing instability for D < 0. Can coupling influence 
this instability? In some works ||^, [58| spinodal shifts 
have indeed been observed. In our model this does not 
occur, for the following reason. By analogy with the feed- 
back mech anism studied above for k = ky, the term in 
Eqn. 

I • 1 our model (unlike 



2.17 that could participate in positive feedback is 
W zz G'(d>)k 2 Sd>. 



§ |38 ) W zz = 



(Eqn. 2.24) so the stability of vorticity fluctuations is 
unaffected by shear. Accordingly, hereafter we consider 
only k = ky. 



SHEAR STARTUP EXPERIMENT 



A. Time-dependence and linear analysis 
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W XZ {t) = Wy Z (t) = 



(5 



(see Fig. |Sp . Although these expressions recover 
Eqn. |2.24| as t ~ + 00 (so that the total shear stress would 
then be on the intrinsic constitutive curve), we show be- 
low that in general the flow goes unstable before this limit 
is reached. 




FIG. 8: Homogeneous background micellar strain W xy vs. t 
for 7 = 2.0, 4.0, 8.0, 16.0 (top to bottom at right of plot). 



The stability analysis of startup flow is more involved, 
because here fluctuations emerge against a background 
state that itself evolves, deterministically, in time. We 
first outline these deterministic kinetics (for an idealised 
noiseless system) before analysing fluctuations. 



2. Inhomogeneous fluctuations 

In a real system, these homogeneous transients repre- 
sent only a background state u(t) = py, W(t), $\, which 
is subject to fluctuations: 



Deterministic "background" kinetics 



At time t = 0, the rheometer plate at y — L is set in 
motion with velocity jLx, giving an instantaneous shear 
rate profile 7(2/, 0) — r yS(y — L). Without noise, the ulti- 
mate steady state would be homogeneous. Firstly, on the 
Reynolds time scale tr = pL 2 /rj the shear rate rapidly 
homogenizes across the cell such that 7(2/) = 7. Secondly, 
on the Maxwell time scale t > tr the micellar strain 
starts to evolve homogeneously, according to Eqn. [I20| 
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(5.2) 

To investigate the stability of these fluctuations, 
we l inear i ze t he model's dynamical equations 
(flf JU7|, (|2~20| ) to get 



d t Su k (t) = M M).6u k (t). 



(5.3) 



This is the counterpart in startup of Eqn. 4.3, with the 
important new feature that M Jt) is time dependent, 
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via its dependence on the homogeneous background state 
u(t) = py> W(i), 4>] and hence on the evolution of the mi- 
cellar strain W_ towards the intrinsic constitutive curve. 
The eigenmodes are therefore now time-dependent: 



W4,a(*)S fe a(*)=M fc (*)]£fc > a(*)- 



(5.4) 



In any startup experiment, then, the micellar strain 
evolves over a time r ss = O(t) (thus defined) towards the 
intrinsic constitutive curve, as described above. The dis- 
persion relation to^^a ( t) co rrespondingly evolves towards 
the one given by Eqn. 4^4 for an initial condition on that 
flow curve. So for a shear rate in the unstable region, at 
least one dispersion branch must go positive at some time 

< T ss so that the homogeneous transient py, W(t). d>] 
goes unstable. In most regimes we find only one posi- 
tive branch | ]6l| and drop the "mode" subscript a, with 
the understanding that we mean the largest branch. At 
wavevector k, the amplitude of the growing fluctuations 
at a time t > to is approximately set by 



A(k, t) ~ exp 



dt'uk(t') 



t 



(5.5) 



We choose a rough criterion for detectability by scat- 
tering to be log A = O(10). This defines a wavevector- 
dependent time scale r lns t(fc), via 



t(fc) 



dt'wkit') = 0(10). 



(5.6) 



Co 



In most regimes, there is a selected wavevector k* at 
which fluctuations emerge fastest, as the result of a peak 
in the dispersion relation £%(t) vs. k. In practice, the 
peak shifts along the k axis in time, but it is still usually 
possible to obtain a reasonable estimate of the overall 
k*; we justify this claim below. We therefore define the 
overall time scale of instability to be 



t(£*). 



(5.7) 



By the time T lnB %, then, the system is measurably in- 
homogeneous, and our linear calculation terminates. In 
general, this occurs well before the intrinsic constitutive 
curve would have been attained, i.e. r inst < t ss (Fig. ||), 
so that the instability is determined not by the time- 
independent dispersion relations of Sec. IV above, but 
by their time-dependent counterparts (given below). 

Is the unstable intrinsic constitutive curve ever at- 
tained before the instability occurs, such that r ss <§; ri nst ? 
A necessary condition is that the growth rate cjf c = 
t^k*(t — t ss ) that would occur once the flow curve were 
reached (given by the dispersion relations of section IV) 
obeys 



Wf c T ss < 1. (5.8) 

This is not usually satisfied (recall Figs. || and ||) since 
u){ c is itself set by the Maxwell time r (with a prefactor 
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FIG. 9: Cartoon: homogeneous startup flow going unstable at 
time rinst before it can reach the intrinsic unstable constitutive 
curve at time r ss . 



set by the slope of the flow curve an d by concentration 
coupling). Nonetheless, condition 5.8 is satisfied just in- 
side the spinodal, since uj[ c — > smoothly at the spin- 
odal. However, this condition is not always sufficient. 



In particular, for shear rates just inside the upper spin- 
odal, the homogeneous micellar strain oscillates strongly 
in startup. Correspondingly, the growth rate significantly 
overshoots u/f c (Figs. |l^ to [l^ below) and fluctuations 
still emerge before the intrinsic constitutive curve would 
be attained. In fact, these oscillations mean that fluctua- 
tions can become (temporarily) unstable in startup, even 
for shear rates above the upper spinodal (as defined via 
slow shear rate sweeps). This upper spinodal is therefore 
not particularly relevant to startup flows. 

For she ar rates just inside the lower spinodal, con- 
dition 5.8 is necessary and sufficient, and the intrinsic 



constitutive curve is then attained before the instability 
develops appreciably. Here we can assume that the sta- 
bility matrix changes discontinuously at t = from the 
stable M (t = 0) (with 7 = F£ = 0), to the unstable 

matrix M (t — 00) for a state on the intrinsic consti- 
tutive curve. The instability is then, even in startup, 
determined by the time-independent dispersion relations 
of Sec. IFvL 



We pause to compare our analysis with that of Cahn 
and Hilliard for a two component system temperature- 
quenched at time t — into the unstable region, 
dfi(4>,T)/d(j) < 0. A good approximation, invariably 
made, is that n(<j>) changes discontinuously at t = from 
its initial stable state to the final one of negative slope, 
i.e. that the heat diffuses out instantaneously with re- 
spect to the time scale at which fluctuations grow. We 
have just seen that the corresponding assumption for 
our purposes (the background state u(t) instantaneously 
reaching the intrinsic constitutive curve) is not in general 
valid. 

In the next section we present results for the time- 
dependent unstable dispersion branch over the time in- 
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terval to — > Tj ns t for several startup quenches, indicating 
in each case the selected wavevector k*. We also give 
results for the time dependent eigenvector (at k*) noting 
whether separation occurs predominantly in the mechan- 
ical variables or in concentration. 



B. Results: uncoupled model 

Fig. [l(^ (top) shows the numerically calculated startup 
dispersion relation Wfe(t) in this uncoupled limit £ — > oo. 
The oscillations arise from the oscillations in Wit) to- 
wards the intrinsic constitutive curve (Fig. |^). Despite 
the time dependence, the main features of the time- 
independent dispersion relation fo r fluc tuations about 
the intrinsic constitutive curve (Fig. 4(a)| ) are still appar- 
ent: there is a Reynolds regime as k — > 0, a non-conserved 
plateau regime at intermediate k, and interfacial cutoff 
at large k. As before, then, in this uncoupled limit there 
is no selected wavevector k* . 

For each startup, we estimated the time Ti nst at which 
the insta bilit y would become measurable, as governed by 
criterion |5.6| applied to wavevectors in the plateau regime. 
It is ma rked by t he thic k line in Fig. [10(a) and an arrow 
in Figs. |1 0(d)| to [l0(f)| . For each value of 7 in Fig. [l(], 
we find Ti ns t <C t ss : instability occurs long before the 
underlying flow curve would have been attained. 

Fig. [l(] (bottom), shows the time-dependent eigenvec- 
tor at wavevector k* — -k (chosen arbitrarily since the 
eigenvector is independent of k in the plateau regime). 
This is dominated by Sj, since SW xy + -qSj = in this 
zero-Reynolds regime. Note also that the normal stress, 
encoded in SZ, dominates the shear con tributi on 5W xy : 
consistently with the remarks of section IV B 1 , the nor- 
mal stress plays an important role in this mechanical in- 
stability. 

The discontinuity in the first derivative of eigenvec- 
tor is due to a crossing of two positive eigenvalues: in 
contrast to the time-independent dispersion relations for 
fluctuations about the intrinsic constitutive curve, in 
startup there is sometimes more than one positive dis- 
persion branch, the smaller one of which can be in the 



subspace [ikSv z ,W xz ,W yz ]. However, this second unsta- 
ble mode only occurs at high shear rates 7 ~ 10, and 
even then only crosses the first for times well after Ti nst : 
consistent with the claim made above, we never observe 
mode crossing in the relevant time regime t < Ti nst . This 
also applies to the coupled model, to which we now turn. 



C. Results: coupled model 

We now give startup results for the coupled model. 
Denoting the experimental DLS value of the diffusion 
coefficient D (table |J) by -D C xpt, figures O and |l2] are for 
D = D oxp t (type I system) and D = lo"D e x P t (type II 
system) respectively. The overall features of these disper- 
sion relations are the same as for their time-independent 



counterparts (Figs. [| and |6(b)| ). In particular there is, 
at any time, a well defined peak k = k pcs ^(t). This peak 
shifts along the k axis in time. At t = to, when co* = 
by definition, we numerically observe that fc pea k = 0. 
However fc poa k very quickly attains a value k* that is 
(practically) time-independent and well approximated by 
Eqn. 4.24 . In this way, the time-dependence of /c pea k 



only occurs at early times t ~ to, for which the growth 
rate is insignificantly small. We argue, therefore, that we 
can choose the ultimate /c poa k = k* as the representative 
wavevector for the instability. 

The time dependent eigenvector at this selected 
wavevector k* is also shown in figures [n] to ||. As 
noted above, the eigenvector encodes the extent to which 
separation occurs in each of different order parameters. 
Experimentally, polarised light scattering is sensitive to 
fluctuations in the micellar strain, while unpolariscd light 
scattering measures fluctuations in the overall micellar 
concentration. In a forthcoming paper p2] | , we explicitly 
calculate the unstable startup static structure factor sep- 
arately for polarised and unpolarised light scattering. In 
this paper, we restrict ourselves to the overall features on 
the instability that are deducible from the eigenvector. 

Fo r type I syst ems at all (unstable) shear rates (fig- 
ures |ll(d)| , ll(c)| an( i H(f)D ) and for type II sy stems 
at shear rates that are not too small (figures 12(e) 



and |12(f)| ), the eigenvector is dominated by the flow vari- 
ables (57 and S W as expected. In contrast, for the type 
II system at low shear rates (Fig. 12(d)) the eigenvec- 
tor is dominated by concentration: here the instability 
is essentially CH demixing, triggered by flow. At higher 
shear rates, even in this type II system, the instability 
is basically mechanical. Note finally that concentration 
coupling affects the relative contributions of the shear 
(W xy ) and normal iZ) stresses. Recall that in the un- 
coupled limit, 5Z ^> SW xy . In this coupled case, the 
shear component can be comparable to SZ (for moder- 
ate applied shear rate 7) or even larger than SZ (low 



VI. CONCLUSION 

In this paper, we have studied the early-time kinetics 
of the shear banding instability in startup flows. Moti- 
vated by recent rheo-optical experiments J3l| in which en- 
hanced concentration fluctuations were observed in the 
flow/flow-gradient and flow/vorticity planes at the on- 
set of instability, we performed a linear stability analysis 
for coupled fluctuations in shear rate, viscoelastic stress 
and concentration using the non-local Johnson-Segalman 
model and a two-fluid approach to concentration fluctu- 
ations. 

We considered two flow histories. The first assumed 
an initial homogeneous state on the intrinsic constitutive 
curve. Using this, we defined the spinodal boundaries 
of the unstable region for slow shear rate sweeps. (Any 
real system can in practice shear band via metastable 
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(a) 7 = 7.0, dispersion relation. (b) 7 = 40.0, dispersion relation. (c) 7 = 70.0, dispersion relation 
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(d) 7 = 7.0, eigenvector at 
log(fc*/7r) = 0.0. 



(e) 7 = 40.0, eigenvector at 
log(fc*/V) = 0.0. 



(f) 7 = 70.0, eigenvector at 
log(fc*/7r) = 0.0. 



FIG. 10: Type A instabilities in a type I system: time dependent dispersion relation (top) and eigenvector (bottom) in the 
uncoupled limit C — * 00 for = 0.11. The rheological model parameters all assume the experimental values of table |. The thick 
line in Fig. 10(a) and the arrows in Figs. d,e,f denote the time at which the instability becomes measurable. The discontinuities 
in the first derivative of the eigenvector components result from a crossing of eigenvalues, discussed in the text. 



kinetics before the unstable region is reached. The anal- 
ogous ambiguity occurs in the CH calculation for fluid- 
fluid, demixing.) For the uncoupled limit of infinite drag 
between micelles and solvent at finite collective micel- 
lar diffusion constant D, fluctuations in the mechanical 
variables (shear rate and stress) decouple from concen- 
tration and are unstable when the intrinsic constitutive 
curve has negative slope, as expected; the concentration 
has its own CH instability for D < 0. 

For finite drag, fluctuations in the mechanical variables 
couple to those in concentration via the positive feedback 
mechanism of Helfand and Fredrickson, and the unstable 
region of what was the purely mechanical instability is 
broadened. In rapid upward stress sweep experiments, 
therefore, "top" jumping should in fact occur before the 
maxiumum in the intrinsic constitutive curve is reached. 
For given values of the coupling parameters, the degree 
of broadening increases with proximity to an underlying 
(zero shear) CH demixing instability. Accordingly, we 
classified systems into two types. Type I systems are far 
from a CH instability, and the mechanical instability is 



only slightly perturbed by concentration coupling. Type 
II systems are close to a CH instability. Type I systems, 
and type II systems at high shear rates, show instabili- 
ties that are predominantly mechanical (type A) . Type II 
systems at low shear rates should show a CH instability 
perturbed by coupling to flow (type B). 

We then discussed the dispersion relations for fluctua- 
tions about the unstable intrinsic constitutive curve. In 
the uncoupled limit, there is a broad plateau in the dis- 
persion relation, with no selected length scale. Concen- 
tration coupling enhances the instability at short wave- 
lengths thereby selecting a wavelength. However the typ- 
ical growth rates predicted by these dispersion relations 
are larger than the rate at which the system can real- 
istically be prepared on the intrinsic constitutive curve. 
We therefore explicitly considered shear startup quenches 
into the unstable region. We assumed that the startup 
flow can be decomposed into a homogeneous background 
state, evolving towards the intrinsic constitutive curve, 
with small fluctuations that (consistently with the above 
remarks) in general go unstable before the the intrinsic 
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FIG. 11: Type A instabilities in a type I system: time dependent dispersion relation (top) and eigenvector (bottom) for a 
coupled model in which all parameters assume the experimental values of table [} The concentration cj) = 0.11. The arrows in 
Figs. e,f show the time at which the instability first becomes measureable (the instability occurs beyond the time window of 
Fig. d). The discontinuities in the first derivative of the eigenvector components result from a crossing of eigenvalues, discussed 
in the text. The instability time Ti ns t occurs beyond the displayed time window for 7 = 2.0. 



constitutive curve can be attained. The main features 
of the time independent dispersion relations for fluctua- 
tions about the intrinsic constitutive curve are nonethe- 
less preserved: the mechanical instability shows length 
scale selection only when coupled to concentration. Our 
results for the time dependent eigenvector at this selected 
length scale confirmed our classification of instabilities 
into types A and B. 

In the coupled model, for small values of the diffu- 
sion coefficient a second lobe of instability develops at 
high shear rates. This could clearly have dramatic con- 
sequences for any putative coexistence of low shear and 
high shear bands, since the high shear band could it- 
self be unstable. Indeed, the high shear band is often 
seen to fluctuate strongly [[33], or to break into smaller 



bands [p3[ . However, in our model this high-shear insta- 
bility is highly sensitive to choice of model parameters, 
and could be an unrealistic feature. 

Our study was confined to fluctuations in the flow gra- 
dient direction, and to the qualitative features of the in- 
stability that can be gleaned from the time-dependent 
dispersion relations and eigenvectors. In a future paper, 
we will present results for the time-dependent unstable 
static structure factor in startup, for fluctuations in the 
entire flow/flow-gradient plane |62| . 
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